ERROR ANALYSIS OF THE BERGMAN KERNEL 
METHOD WITH SINGULAR BASIS FUNCTIONS 



M. A. LYTRIDES AND N. S. STYLIANOPOULOS 

Abstract. Let G be a bounded Jordan domain in the complex plane 
with piecewise analytic boundary. We present theoretical estimates and 
numerical evidence for certain phenomena, regarding the application of 
the Bergman kernel method with algebraic and pole singular basis func- 
tions, for approximating the conformal mapping of G onto the normal- 
ized disk. In this way, we complete the task of providing full theoretical 
justification of this method. 



1. Introduction 

Let G be a bounded, simply-connected domain in the complex plane C 
whose boundary T := dG is a Jordan curve and let O := C \ G denote the 
complement of G with respect to the extended complex plane. Fix zq € G 
and let /o denote the conformal map of G onto the disk D(0, tq) := {z : 
\z\ < ro}, normalized by the conditions fo{zo) = and /o(2o) = 1- The 
quantity ro := ro{G, zq) is called the conformal radius of G with respect to 

Zq. 

For the inner product 



{f,g) := / fiz)giz)dA{z), (1.1) 
Jg 

where dA denotes the differential of the area measure on C, we consider the 
Hilbert space 

Ll{G) :={/:/ analytic in G, (/, /> < oo}, (1.2) 

with corresponding norm ||/||l2(g') := (/, /)2. 

Let K{-,zq) denote the Bergman kernel function of G with respect to zq. 
This is the unique function of L'^{G) satisfying the reproducing property 

{g, K{; zo)) = g{zo), for ah g G LI{G). (1.3) 
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2 M. LYTRIDES AND N. STYLIANOPOULOS 

It follows from (ll.3|] that the kernel K{-,zo) is related to the mapping func- 
tion /o by means of 

^«<'> = l(fe^ '"^ AU) = ^^/V(C.^«W. (1.4) 

see e.g. [3l p. 33]. These yield the two relations, 



K{z,zo) = ^foiz) and rp = ^^^ =. (1.5) 
vrr^ V7rK(zo,2o) 

Now, let {P„}J^o denote the sequence of the Bergman polynomials of G. 
This is defined as the sequence of polynomials 

Pn{z) = XuZ"" + ■ ■ ■ , Xn>0, n = 0,l,2,..., (1.6) 

that are orthonormal with respect to the inner product (jl.ip . i.e., 

I Pm{z)Pj^dA{z) = 5,n,n. (1-7) 
JG 

The Bergman polynomials form a complete orthonormal system in L'^{G). 
Therefore, in view of the reproducing property ()1.3p . 



K{z,zo) = yPj{zo)P,{z), {U 



j=0 



locally uniformly with respect to z £ G. 

The Bergman kernel method (BKM) is an orthonormalization method for 
computing approximations to the conformal map fo{z). It is based on the 
fact that the kernel K{z,zq) is given explicitly in terms of the Bergman poly- 
nomials {Pn{z)}^^Q. Thus, the partial sums of the Fourier series expansion 
of K{z, zq) are given by 

n 

Kn{z,zo) ■.= J2W^Pii^)^ n^^- (1-9) 

j=0 

The polynomials {Kn{z, zo)}^^q are the so-called kernel polynomials of G, 
with respect to zq. They provide the best L?'{G)- approximation to K{-, zq) 
out of the space P^i of complex polynomials of degree at most n. 

In accordance with (jl.4p . the n-th BKM approximation to /o is given by 

^n{z):=- \ r/ K„„i(C,2o)f^C, n^n. (1.10) 



Kn-l{zQ,ZQ) 

This defines the sequence {tt^}^^ of the Bieberbach polynomials of G, with 
respect to zq. The polynomial 7r„ solves the following minimal problem: Let 

K-={P-P^ P„, with p{zo) = and p'(zo) = !}■ 
Then, for each n G N, the polynomial Tin minimizes uniquely the two norms 
ll/o -p'\\l^{G) and ||p'||l2(g) over all p G P^; see e.g. O Kap. Ill, §1]. 

Regarding the convergence of the method, we note that in cases when /o 
has an analytic continuation across T, then this is a consequence of Walsh's 
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theory of maximal convergence |22^ §4.7, §5.3]. In order to be more specific, 
let ^ denote the conformal map of O onto A := {w : \w\ > 1}, normalized 
so that near infinity, 

=7z + 7o + — + ^ + --- , 7>0. (1.11) 

Note that 7 = l/cap(r), where cap(r) denotes the (logarithmic) capacity 
of r. Then, 

ll/0-^n|Loo(G)=o(;^), (1.12) 

holds for any 1 < i? < |<I>(zi)|, but for no R > |<I>(zi)|, where zi denotes 
the nearest singularity of /o in $7; see also [3l Ch. I]. (We use || • H^ooj-fj) to 

denote the sup- norm on G.) 

In cases when T is piecewise analytic and /o has singularities on F, then 
Levin, Papamichael and Siderides were the first to observe in jTl] that the 
error (I1.12P depends on the boundary singularities of the mapping func- 
tion /o on r, and also on the singularities of the extension of /o across the 
segments of T into il.. Accordingly, in order to improve the numerical per- 
formance of the BKM, they extended the method by orthonormalizing a 
system of basis functions consisting from monomials, as in the BKM, and 
also from functions that reflect the dominant singularities of /o on T and in 
Q. This extension is known as BKM/AB (AB stands for augmented basis). 
The BKM/AB was used subsequently in [15] and [16j . 

The most precise results regarding the convergence of the BKM are due 
to D. Gaier [6]. In particular, under the assumption that F is piecewise 
analytic without cusps, Gaier derived the estimate 

ll/o - 7r„||^oo(G) = C'(logn)-^, (1-13) 

where s := A/ (2 — A) and Avr (0 < A < 2) denotes the smallest exterior angle 
where two analytic arcs of T meet. Regarding sharpness of the estimate 
(jl.l3p . it was shown in [U Thm. 4] that there are cases where the exponent 
s can not be replaced by a smaller number. However, the factor logn can 
be replaced by \/log n, see [T] and |12( Rem. 3.1]. A lower estimate of the 
form 

||/o-vr„||i^(G) >c— , (1.14) 

provided that 1/(2 — A) is not a positive integer, where c is a constant that 
does not depend on n, was established in [12^ Thm. 3.2] by Maymeskul, Saff 
and the second author. 

The theoretical justification of the BKM/AB with basis function that 
reflect the corner singularities of /o was given in [12], by means of sharp 
estimates for the associated BKM/AB errors. The purpose of the present 
paper is to derive theoretical results that justify the use of basis functions 
that reflect (a) pole singularities of /o and (b) both corner and pole singu- 
larities of /q. More specifically, we derive upper and lower estimates for the 
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BKM/AB errors in the case (a), and upper estimates for the BKM/AB er- 
rors in the case (b). In doing so, we complete the task that was put forward 
by Yu. E. Khokhlov, reviewer of the introductory paper [TT] of the BKM/AB 
in the Mathematical Reviews, who concluded that: proof of the conver- 
gence of the numerical method given and an investigation of its convergence 
rate are lacking, so the results obtained are of a heuristic nature" . 

The paper is organized as follows: In Section 2 we set up the notation and 
recall the BKM/AB. Section 3 is devoted to the study of the various BKM 
and BKM/AB errors, in cases when /o has an analytic continuation across 
r, hence only basis functions reflecting poles are used in the BKM/AB. 
In Section 4, we consider the case when both corner and pole basis func- 
tions are included in BKM/AB. Finally, in Section 5, we present numerical 
computations that illustrate the theory of Sections 3 and 4. 

2. The Bergman kernel method with singular basis fungtions 

2.1. Corner singularities. Throughout this section we assume that the 
boundary F of G consist of N analytic arcs that meet at corner points r^, 
k = 1,2, . . . , N , where they form interior angles afcvr, < < 2. Then, we 
have the following asymptotic expansions for /q, valid near r^: 

(i) If Ok is irrational, then 

foiz) = /o(rfc) + Bp,g{z - rfcf+''/"^ (2.1) 

where p and q run over all integers p > 0, q > 1 and i?o,i 7^ 0. 

(ii) If Ofc = a/b, with a and b relative prime numbers, then 

fo{z) = foiTk) + Bp,,,m{z - rfc)f+'?/"Hlog(z - Tk)r, (2.2) 

p,q,m 

where p, q and m run over all integers p>0,l<q<a,l<m<p/b 
and ^0,1,0 / 0. 

(iii) If Tfc is formed by two straight-line segments, then 

oo 

fo{z) = foin) + Y Bi{z - nf"", (2.3) 
1=1 

where Bi ^ 0. Furthermore, (|2.ip holds in the case when is 
formed by two circular arcs, or a straight-line and a circular arc. 

In the above, (i) and (ii) are due to Lehman [9], while (iii) emerges easily 
from the reflection principle; see also [U §2.1] and [Ml pp. 6-7]. 

It follows from (iii) that if G is a half-disk or a rectangle, then /o has 
a Taylor series expansion valid around each corner, and thus an analytic 
continuation across F into i}. In this case, the only singularities of /o are 
simple poles in il.. This shows that the study of the BKM/AB, even with 
only pole basis function is important in the applications. 
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For simplicity in the exposition, we shall assume throughout this paper 
that no logarithmic terms occur in the asymptotic expansion of /o near 
the corner r^, k = 1,2, . . . , N . This, for example, will be the case in the 
expansions (I2.ip and (I2.3P above. Nevertheless, our method of study can be 
adjusted to cover logarithmic singularities as well. 

Let M denote the number of corners of T for which ak is not of the special 
form 1 /m, m G N. When we present results for corner singularities we shall 
assume that M > 1. We index such corners by r^, k = 1, . . . , M. That is, if 
N > M, then the mapping function /q has an analytic continuation in some 
neighborhood of the corner Tjy. 

For k = 1,...,M, we denote by {7^^^^}^;^ the increasing arrangement 

of the possible powers p + q/ak of (z — Tk) that appear in the asymptotic 

expansion of /o(-z) near r^. In particular, if Tk is formed by two straight-line 

segments, then 7^^^^ = j/uk, j = 1,2,.... Also, if ak is irrational, or the 

corner Tk is formed by two circular arcs, then 

W 1 / 

7| = 1/ak; 

7^^^ = 1/afe + min(l/afc, 1); 

1/ak + 2, 0<ak< 1/2, 
7^^ = <; 2/ak, l/2<Qfc<l, 
1/ak + 1, 1< a/c < 2; 



Remark 2.1. Under the assumption regarding the no-appearance of loga- 
rithmic terms, the asymptotic expansion near Tk, k = 1,2, .. . ,M, can be 
written in the form 

00 

uz) = Y,4\z-Tkr^\ (2.4) 

j=0 

where, j^^^ := and a^^^ ^ 0. Note that, we always have 7!^^^ > 1/2, and 

(k\ I , > (fc) 

since Tk is not a special corner, 7]^ ^ N. Therefore (z — Tk)"^^ has an 

algebraic singularity at r^. However, when ak is rational, it is possible that 
7]^^^ E N, for indices j > 2, so that (z — Tk)'^^ is analytic at r^. 

2.2. Pole singularities. Since /o(zo) = 0, zq S G, it follows from the 
reflection principle for analytic arcs that the extension of /o across any 
segment constituting F would have a pole or a pole-type singularity at the 
reflected images of zq. For example, if F consists explicitly from straight- 
line segments and/or circular arcs, then /o has a simple pole (due to the 
univalency of /o) at every mirror image of zq (with respect to the straight- 
lines) and at every geometric inverse of zq (with respect to the circular arcs) , 
that lies in $7. More generally, /o may have at points z-,- € a pole, or a 
poly-type, singularity of the form 

(z-ZJ)"^^/"^ fcj,mjGN. (2.5) 
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According to [TTl § 5.1], the following three special cases occur frequently in 
the applications: 

(i) kj = iTij = 1. In this case, /o has a simple pole at Zj. 

(ii) kj = 2, nij = 1. In this case, /o has a double pole at Zj. 

(iii) kj = 1, rrij = 2. In this case, /o has a rational pole singularity at Zj. 

In order to describe the BKM/AB, we assume that the nearest singu- 
larities of /o in $7 are poles or rational poles, of the form ()2.5p at points 
Zj, j = 1,2, ...K, where \^{zi)\ < |$(z2)| < ••• < and that the 

other singularities of /o in Q, occur at points Zk+i,Zk+2, ■ ■ ■, where < 

mzK+i)\ < mz,+2)\ < .... 

2.3. BKM/AB. Using the above notation, the BKM/AB with n mono- 
mials, K poles and corner singularities at each (non-special) corner r^, 
k = 1,2 ... , M, can be summarized as follows: 

(i) Start with the augmented system {r]j} consisting of: 

(1) the nearest poles or rational poles, i.e., for j = 1, 2 . . . k, 



z — z. 



(2.6) 



(2) the dominant tm '■= X^fcliPfc algebraic singular functions, i.e., 
for each non-special corner Tk, k = 1,2, . . . , M, 

VK+jiz) = 

(3) the monomials 



7]^+jiz) = [{z-Tkr'f\', j = l,2,...pk; (2.7) 



VK+rM+ji^) = i^^) ^ i = l,2, ...,n. (2.8) 

(As it was noted in Remark 12.11 it might be possible that 7^*^^ G N. 
If this happens, we avoid redundancy in the basis by omitting such 

7f.) 

(ii) Orthonormalize {rjj}, by means of the Gram-Schmidt process to pro- 
duce the orthonormal set {Pj}, where 

_ j 

P,{z) = Y,hj,^r^i{z), bj,j>0. (2.9) 

1=1 

(iii) Approximate K{z, zq) by its finite Fourier expansion with respect to 

Kn{z,zo):= Pji^o)Pj{z)= dnjrjjiz). (2.10) 

(iv) Approximate fo{z) by 

7f„+i(z) := -/ Kn{C,zo)dC= Y (^njfJ'jiz), (2.11) 

KniZf), Zn) Jzn -I 
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where 

fi,{z) := /%,(C)dC. (2.12) 



We call the functions {Pj} the augmented Bergman polynomials of G, with 
respect to {?7j}, and the functions {vTn} the augmented Bieberbach polynomi- 
als over the system {nj}- Clearly, TTn{zo) = and tt'^{zo) = 1, n G N. Note 
that {Pj}'^i forms a complete orthonormal system in L^(G). Consequently, 



Kiz,zo) = Y.Pjizo)Pjiz), (2.13) 
i=i 

locally uniformly with respect to z £ G, cf. (jl.Sp . 

We conclude this section by presenting a result which shows that the two 
errors H/q — TTn+i II I/2(g) and ll-^('; •2;o) — -fCn(', 2o)||i2(G') are of the same order. 
This fact will be used below in Sections 3 and 4. 

In what follows we denote by c, ci, C2, . . . , constants that are independent 
of n. For quantities A > 0, > 0, we use the notation A ^ B (inequality 
with respect to the order) A < cB. The expression A ^ B means that 
A ^ B and B ^ A simultaneously. 

Lemma 2.1. It holds that, 

ll/o - K+iWlHG) ^ \\K{; zo) - kr,{; zo)||l2(g)- (2.14) 

Proof. We set m:= k + tm + n and note that ([131), (f230D - (l2T3D . imply: 



/^(z)-^;+i(z)=vrr2 5]P,(zo)P,(z)-{5^|P,(zo)lVE^:'(^o)P,(^ 
j=i j=i j=i 



: [vrrg - {J2 \P,{zo)f}-']P,{zo)P,{z) 
i=i i=i 



+7Trl J2 Pj(^o)Pj{z). 

j=m+l 

Therefore, by using the orthonormality of Pj we see that, 

m oo 

m-K+iWhiG) = E [^r^o-'^/Kn{zo,zo)]'\Pj{zo)\'+{^rif ^ |P,(^o)|'. 

j=l j=m+l 

Now, using once more (|1.5|) . we obtain, after some trivial calculation, that 

ll/o - 7rn+llli2(G) = K{zo,Zo) - kn{zo,Zo) K {zq, Zq) kn{zo, Zq) 

This and (jl.3p . with (7(-) = K{-,zo) — Kn{-,ZQ), leads to 

II f ^' l|2 ll^(-,^0)-^n(-,^0)lli2(g) 

ll/o - ^n+lllL2(G) = — ^ (2.15) 
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and the result ()2.14p follows from the set of the obvious inequalities, 

\Pi{zo)\ = Ki{zo,zo) < Kn{zo,zo) < K{zo,zo) = (7rro)"\ 

with constants depending on tq and |Pi(zo)| only. □ 

Remark 2.2. It is clear from the proof that the result of Lemma l2.ll holds 
true for any complete orthonormal system. We note that for the system 
{Pj}'j^i to be complete it suffices that F is a bounded Jordan curve. In 

particular, (I2.14p holds with vr^+i and Kn in the place of TTn+i and Kn. 

3. BKM/AB WITH POLE SINGULARITIES 

In this section we study the BKM and BKM/AB errors ||/o~'''"n|lioo(f;) and 
ll/o ~^nllLoo(G)' u^der the assumption that /o has an analytic continuation 
across T in tl and its only singularities are poles, or rational poles, of the 
type ()2.5p . More precisely, we refine the classical result (jl.l2p for the BKM 
error, and at the same time we obtain a lower estimate for it. Furthermore, 
we establish upper and lower estimates for the BKM/AB error. The lower 
estimates and the refinement are obtained by exploiting the assumption 
regarding the singularities of /o and by using certain important results of 
E.B. Saff on polynomial interpolation of meromorphic functions |I8]. Since 
the results of [18] were established for domains with smooth boundaries, we 
show in the next lemma that they hold true for domains with corners. 

In order to do so, we use the Faher polynomials {Fnj^^Lo ^- recall 
that Fn{z) is defined as the polynomial part of the Laurent series expansion 
of at infinity, i.e., 

- = of i ), as 2 ^ oo. (3.1) 



This, in view of (jl.lip . gives Fn G Pn and 

F„(z) = 7"^" + -- - . (3.2) 

Let L/? {R > 1), denote the level curve of index R of <I>, i.e., 

Lr:={z:\^[z)\=R], (3.3) 

so that Li = F. Note that Lr, for i? > 1, is an analytic Jordan curve. We use 
Gr to denote its interior, i.e., Gr := \ni{LR). The following result gives the 
exact rate of convergence of the minimum uniform error in approximating 
meromorphic functions by polynomials. 

Lemma 3.1. Assume that the boundary F of G is piecewise Dini-smooth 
and consider a function f which is analytic on Gg, for some q > 1, apart 
from a finite number of poles on Lg. Let m denote the highest order of the 
poles of f on Lg. Then, 

m-l 

-/Jl/-Pll.»,g,-— . (3.4) 
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A curve T is piecewise Dini-smooth if it consists of a finite number of Dini- 
smooth arcs. An arc z = z{s), where s £ [a, 6] stands for the arclength, is 
caUed Dini-smooth if z'{s) is continuous on [a, b], and if z'{s) has a modulus 
of continuity uj which satisfies /q [uj{t)/t] dt < oo, for some a > 0. We note, 
in particular, that a piecewise Dini-smooth curve may have corners or cusps 
and that a piecewise analytic Jordan curve is also piecewise Dini-smooth. 

Proof. We recall the following two facts regarding Faber polynomials: 
(i) For any r, R, with 1 < r < i?, it holds 



F„(z) = <i>"(z)|l + 0(^— J|, zeLn, (3.5) 

see e.g. [HI p. 43]. 

(ii) Under the assumption on F, the Faber polynomials are uniformly 
bounded on G (see 0), i-e., 

lli^nllioo(G) <c(r), nGN, (3.6) 

where c(F) is a positive constant that depends on F only. 

Observe that ()3.5p implies that the sequence {Pn{z)}^^i has no limit 
point of zeros exterior to G. Also, from (j3.5p and ()3.6p we have for z £ G 
and t G Lg that, 

'^-(")'<^, nEN. (3.7) 



l^n(t)| - 

Now, following the proof of Theorem 2 of [18] and using the sequence of the 
Faber polynomials {Rn} in the place of {w^}, we conclude that there exist 
polynomials {pn}^=i, such that 

^m— 1 

ll/-PnllL-(G) <C2(r)^^, neN, (3.8) 

see also [19^ p. 399]. This yields the upper bound in (13. 4p . The lower bound 
follows at once from Theorem 10 of p^, by observing that is simply- 
connected and hence its Green function with pole at infinity has no critical 
points. □ 

The following result is the so-called Andrievskii's lemma for polynomials 
and rational polynomials. Its proof, for bounded Jordan domains such that 
the inverse conformal map g : O ^ G satisfies a Lipschitz condition on 
D, can be found in [4]. This condition is certainly satisfied by the type of 
domains considered below. 

Lemma 3.2. Assume that F is piecewise analytic without cusps. Then: 
(i) Por any Pn G Pn, with Pn{zo) = 0, it holds 

ll^'nllLoo(G) ^ \/logn||P^||i2(G), n > 2. (3.9) 
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(ii) For any Pn G Fn, with Pn{zo) = 0, and q a fixed polynomial with no 
zeros on G, it holds that 



WPn/qh^^G) ^ V^WiPn/qyU^iG), n > 2. (3.10) 

3.1. BKM. The next theorem complements the classical result (jl.l2p of 
Walsh, in the sense that it provides a lower estimate and, in addition, uses 
the precise g = \^{zi)\ in the denominator, instead of any R, with 1 < R < g. 
This is done by utilizing extra information on the nature of the singularities 
of /o in 0. 

Theorem 3.1. Assume that T is piecewise analytic without cusps. Assume 
further that the conformal map /o has an analytic continuation across T, 
such that /o is analytic on Gg, for some g > 1, apart from a finite number 
of poles on Lg. Let m denote the highest order of the poles of /o on Lg. 
Then, 

n^~^ \/log n 

—;^^\\f0-TTn\\L^(G)^ ' ^^^2. (3.11) 

Proof We observe first that the kernel K{z,zq) shares the same analytic 
properties with /o on Gg, apart from an unit increase on the order of its 
poles on Lg. Therefore, using Lemma [3.11 with / = K{-,zo), we conclude 
that 

\\K{-,Zo)-pn\\L^(G)^—' (3-12) 

for some sequence of polynomials {pn}^=i- Since the L^(G)-norm is domi- 
nated by the L°°(G)-norm, (I3.12p leads to the estimate 

\\K{-,Zo) -Pn\\L2{G) ^ — • (3-13) 
Then, the minimum property of the kernel polynomials implies that 

\\K{;Zo)-Kn{;Zo)\\mG) ^ (3.14) 

Q 

which, in conjunction with Remark 1 2. 2 [ yields the estimate 

\\fo-<\\LHG)^^- (3.15) 

Now, we use Andrievskii's Lemma l3.2f i) and employ the method of An- 
drievskii and Simonenko, see e.g. O §2.1]. This method enables the transi- 
tion from an upper bound of the error ||/o — 7r^||L2(G) ^ similar bound for 
the error ||/o ~7i"n||x,c3o((j), with the extra cost of a \/log n factor, and leads to 
the upper estimate in ()3.1ip . The lower estimate follows immediately from 
Lemma 13. 1[ □ 



The following pointwise estimate is useful in the study of the distribution 
of the zeros of the Bergman polynomials; see e.g., [10], [H], [19] and [8]. 
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Corollary 3.1. With the assumptions of Theorem \3.1\ it holds, 

\Pn{zo)\ n G N, (zo G G). (3.16) 

Proof. The result emerges easily from (|3.14p , using the reproducing property 
of K{-^zo), the fact that Pn+i is orthogonal to any polynomial 
the Cauchy-Schwarz inequality: 

|P„+l(zo)| = |(P„+l,/^(-,Zo))| = K^n+l,if(-,2o)-if„(-,Zo)>| 
< ll^n+l||L2{G) \\K{-,Zq) - Kn{-,ZQ)\\L2(^G) 
= \\K{-:Zq) - Kn(-,Zo)||L2(G). 

□ 

3.2. BKM/AB with pole singularities. We exploit now the specific as- 
sumptions on the singularities of the analytic extension of /o studied in 
Section 12.21 More precisely, the assumption that the nearest singulari- 
ties of /o are k poles, each one of order kj at Zj, j = 1,2, ...k, where 
|^(^i)| < l'^'(-22)| < ••• < I*^(-2k)|, and that the other singularities of /o 
occur at points z«;+i, 2;k+2, • • •, where |^>(2;k)| < |$(2;k+i)| < \^{zi^+2)\ ^ • • • • 
Therefore, for the BKM/AB we consider the system defined by the sin- 
gular functions in (j2.6p . with rrij = 1, j = 1, 2, . . . , k, and the n monomials 
in (j2.8p . Accordingly, we let P^^ denote the following space of augmented 
polynomials: 

K+n 

Pn^ := {p : p{z) = t,rij{z), t, G C}. (3.17) 

We note that the associated augmented kernel polynomial Kn{z,ZQ) is the 
best approximation to K{z,zq) in L'^{G) out of the space P^^, i.e., 

\\K{-,zo) - Kn{-,zo)\\L2(^G) < \\K{-,zo) -p\\l^g)i (3-18) 
for any p G P^^ . 

The next theorem provides an estimate for the error in the resulting 
BKM/AB approximation to /q. 

Theorem 3.2. Assume that T is piecewise analytic without cusps and set 
g := |$(zk+i)I- Then, 

Wfo-^nWi^fG)^^, (3-19) 
for any R, with 1 < R < g, but for no R > g. 

Proof. Observe that K{-,zo) has poles of order kj+1 at each Zj, j = 1,2, ... , k, 
and set Q{z) := nj=i('2 ~ Zj)^^~^^. Then, the function K{z, zq)Q{z) is ana- 
lytic in the interior Gg of the level curve of Lg, and from Walsh's maximal 
convergence theorem [22^ §4.7] it follows that, for any R, with 1 < R < g, 
there exists a sequence of polynomial {pnj^i, such that 

\\K{-,zo)Q-Pn\\Loo^G)^^^ ™eN. (3.20) 
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Let now d := mmj=i^2,...,K{k ~ '■ z £ T} denote the distance of F from 
the poles {zj}^^^, and set ^ := X]j=i%- Then, |Q(z)| > d'^'^^, z G T, and 
(IM]) gives 

WTc-f \ ^"11 _ ^ C 1 

Since the L^(G)-norm is dominated by the L°°(G)-norm, we see that there 
exist a sequence of rational polynomials {QnlJ^Ln with Qn G F^^, such that, 

\\K{-,zo) -Qn\\L2(G) ^ nGN. 

Therefore, using the minimum property (j3.18p of the augmented kernel poly- 
nomials, we have 

1 



\\K{-,zo)-Kni-,zo)\\LHG)^j^, ^eN, (3.21) 
and this, in conjunction with the equivalence Lemma [2. 11 yields the estimate 



ll/o - ^nllL2(G) 



1 



n G N. 



Next, we recall that 



-n-niz) 



{z — ZjY^ 



{zq 



Zjfj 



I.e., 



vr„(2;) 



where ^(z) := Y\_i^ ~ ^i)^^ 



(3.22) 



(3.23) 



(3.24) 



i=i 



and ^(2;) is a polynomial of degree n + 1^. 

Then, the transition from the L^(G)-norm to the L'^(G)-norm is done as 
in the proof of Theorem 13. H where now, in view of ()3.24p , Lemma I3.2( ii) is 
applicable. This leads to, 

VTogn 



ll/o 



■^"IIl°°(G) ^ 



n> 2, 



(3.25) 



and (j3.19p follows with a different, but still less than g, R. 

Finally, the fact that (j3.19p holds for no R > gis evident from |22l Thm. 6, 
Ch. IV], since the contrary assumption would lead to the contradicting con- 
clusion that /o has no singularities on L^; see the next remark. □ 

Remark 3.1. From (j3.23p it is clear that TTniz) = qniz) +Pniz), where is 
defined by the nearest k poles of /o in Q and pn G Pn- Hence, (j3.19p gives 

1 



IL°°(G) - 



II (/o - ^k) -Pr, 

for any 1 < R < g, and Theorem 6 in |22l Ch. V] implies that the function 
/o ~ Qk is analytic in Gg. This shows that the the rational polynomial q^, 
constructed by the BKM/AB considered above, cancels out the specific poles 
of /o that contains. In particular, this provides the theoretical justification 
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for the heuristic observation made to that effect by Papamichael and Warby 
in [m p. 652]. 

A finer estimate than ()3.19p can be obtained if the singularities of /o on 
-^|$(2„+i)| are a finite number of poles. 

Theorem 3.3. Assume that T is piecewise analytic without cusps and set 
g := |<I>(2:k+i)| ■ Assume, in addition to Theorem \3.2l that /o has a finite 
number of poles and no other singularities on Lg and let m denote their 
highest order. Then, 

— ^ ^ ||/o-vr„||i«.(G) ^ , n>2. (3.26) 

Proof. The upper estimate follows by working in the same way as in the 
proof of Theorem 13.21 but using here the precise result of Lemma [3.11 in the 
place of Walsh's theorem in (j3.20p . 

To obtain the lower estimate, observe that qTTn is a polynomial of degree 
n + ^ (see (|3.24|) ) and that the function qfo is analytic on G^, apart from a 
finite number of poles on Lg. Hence, from Lemma l3. II we have, for n G N, 

Ik/o - g^nllL^(G) ^ inf \\qfQ-p\\^^ -^ y——, (3.27) 

which yields the estimate 

\\fo-^n\\L^(G) > _ 

II9IIl°°(g) S 

and hence the required result. □ 

In the more general case, where the nearest k singularities of /o in il. are 
rational poles of the type ()2.5p . we have the following result regarding the 
associated kernel polynomials Kn{-,zo). 

Theorem 3.4. Assume that T is piecewise analytic without cusps and set 
g := |$(zk+i)|- Then, 

1 

for any R, 1 < R < g. 

Proof Set Q{z) := U'j=i{z- Zjp/"^^+^ and follow the proof of Theorem EJ 
up to Equation ()3.2ip . □ 

4. BKM WITH POLE AND CORNER SINGULARITIES 

In this section we assume that /o has a singularity on T and study the 
BKM and BKM/AB errors, corresponding to a variety of different syntheses 
of the system {rjj} of basis functions. In stating the results we use the 
notation and the assumptions set up in Sections 12.11 and 12. 2[ 



^m-l 



\K{;Zo)-Kn{;Zo)\\L2iG) ^ 1^, (3-28) 
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4.1. BKM. Our first result is a straightforward consequence of Theorem 3.1 
of |19j and Lemma 13.21 above. 

Theorem 4.1. Assume that T is piecewise analytic without cusps and set 
g := |$(zi)| and s := min{(2 — OLk)/ak '■ ^ <k < A/}. Then, 

ll/o - 7r„||^,x,(G) < ci\/iog^;^ +C2-^, n>2. (4.1) 
for any R, I < R < g. 

Proof. Observe that in Theorem 3.1 of [19] can be chosen arbitrarily close 
to Zi. Thus, from the minimum property of the kernel polynomials Kn{-, zq) 
we have 

\\K{-,zo)-Kni-,zo)\\L2(G)<ci^ + C2^, n E N, (4.2) 

for any R, 1 < R < g, and the transition from the L^(G)-error in ()4.2p to the 
L°°(G)-error in (14. ip . goes in the same lines as in the proof of Theorem 13.21 

□ 

Remark 4.1. Clearly, as re — )• oo, ()4.ip yields the result ()1.13p . However, 
Theorem 14. 1 1 does more: It captures, in a very precise form, the dependance 
of the BKM error ||/o — ^n||£,oo(g) for "small" values of re, on both the corner 
and pole singularities of /q. This dependance has been testified numerically 
in [11] and has given rise to the introduction of the BKM/AB. 

The following result is a simple consequence of (j4.2p . Its proof is similar 
to that of Corollary 13.11 

Corollary 4.1. With the assumptions of Theorem \4.1\ it holds, 

\Pn{zo)\ < ci^ + C2^, n G N, {zoG G). (4.3) 
re^ R^ 

Remark 4.2. Since |P„(2:o)| < ||i(r(-,2:o) — Kn{- , zq)\\ ^2 (^q-^ , it follows from 
Corollary 14.11 that, if for small values of re, Pn{zo) decays geometrically to 
zero, then the most "serious" singularity of K{-, zq), and hence of /o, is the 
nearest pole in and not an algebraic singularity on the boundary, as the 
asymptotic estimate (|1.13p would suggest. On the other hand, given that 
/o has a singularity on T, Theorem 2.1 of [1^ implies that any point of T 
is a point of accumulation of the zeros of the sequence {Pn}^=i- Therefore, 
an easy way to check whether a pole singularity is more serious than an 
algebraic singularity, for a range of values of re, is by plotting the zeros of 
Pn for the same range: // the zeros stay away from a specific part of the 
boundary, this indicates that Pnizo) decays geometrically and therefore the 
presence of a pole singularity near that part. We refer to [191 Examples 2, 
3], where (14. 3p was used as the tool for explaining the misleading nature of 
such plots. 
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4.2. BKM/AB with corner singularities. From our assumptions on F, 
it follows that the conformal map /o can be extended analytically, by means 
of the reflection principle, beyond F to a larger Jordan domain G, such that 
the boundary dG of G consists of analytic arcs to be fixed below. For this, we 
recall our assumptions on the position of the nearest poles Zj, j = 1, . . . , k, 
of /o in $7 and pick up a point Ci near zi, but interior to the level curve 
Lg, with Q := \^{zi)\. Next, we draw the level curve Lg, with g := |$(Ci)| 
and fix on it points (k, k = 2, . . . , N , "between" and Tk+i, where we set 
ttv+i = Ti. We connect each non-special corner r^, fc = 1, . . . , M, with the 
two adjacent CkS, by using two analytic arc. Next, we denote by l/. the two 
arcs emanating from and call Z^r the part (or parts) of the level line 
that joins together those consecutive points Cfc that have only one connection 
with Tk- See Figure [H for a possible arrangement of corners r^, points Cfc; 
and arcs 1^ and lj\f. Finally, we define G by taking dG := {^^Lih} U 1^. 
The above construction is such that: 

(i) dG is a piecewise analytic Jordan curve that meets F at the non- 
special corner r^, fc = 1, . . . , M. 

(ii) /o is continuous on GUdG and analytic in G and on dG, except for 
the endpoints r^. 

(iii) The asymptotic expansion (j2.4p holds for z G l^, k = 1, . . . ,M, in 
the sense that, for any Pfc € No, 

Pk 

/o(z) = X;af (z-r,.)^r+7(.) Jz), 

j=o (4.4) 

7w ^ (z) = of(z-Tfc)^^'+^ 

We consider now the application of BKM/AB with only corner singulari- 
ties, where we use pk S Nq singular function for each non-special corner r^, 
k = 1,2, . . . , M. In order to measure the BKM/AB error we set 

Uk := min{i > pk : 7f ^ ^ N, af^ / 0}, (4.5) 

and assume that at least one of z/^'s is finite, otherwise the results be- 
come trivial. The associated BKM/AB system {rjj} is thus defined by 
'"Af = X^fcLiPfc singular functions of the form (\2.7h and n monomials (\2.8h . 
Accordingly, we let P^^ denote the space of augmented polynomials: 

Clearly, the associated augmented polynomial Kn{z,zo) is the best approx- 
imation to K{z, zo) in L'^{G) out of the space P^^. 

Let TTn denote the BKM / AB approximation resulting from . Then we 
have the following: 




Theorem 4.2. Assume that F is piecewise analytic without cusps and set 
Q := |$(zi)| and s* := min{(2 - ak)-fi''J : I < k < M}. Then, 

ll/o -7rn|lioo(G) < ci\/logn-^ + C2-^, n > 2, (4.7) 
for any R, 1 < R < q. 

Proof. Using Cauchy's integral formula for the derivative of the extension 
of /o we have, for z & G, 



m = ^.Lj!^,dt 

2vrz Jqg [t- zy 

For each Tk, k = 1, . . . ,M, we consider the first terms up to pk, of the 
Lehman expansion (j4.4p for /q: 

(fc) 

Ffe(z):= J^af (z-Tfc)^.- . (4.9) 
i=o 
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Since the function Fk{z), is analytic in G and continuous on dG we have, as 
in gS]), for z£G, 



Therefore, 



■ " - , 1 /■ F,(t) 



fc=l fc=l r=l ^ ^ fc=l -^'^ 



/, _ y\2 Ott^- 



27ri ^ /,, (t - z)2 27ri ^ ^ J, it- z)' 

1 ^ r Fk{t) , 



■dt 



Hence, for z G G, 

M 



f[,{z)-^Fl{z)=g{z)+h{z), (4.10) 



k=l 



where, 



and 



27riJ,{t-z)^-^" 2m^^^Ji^{t-zY 



-E 



(4.12) 



^ (t - z) 



Now, we denote by r = 1, . . . , M, the part of the level hne Lgr that 
shares the same endpoints with Ir, so that = {U^il[.} U Itv and Ij. U 1^ is 
the boundary of a Jordan domain in see Figured! Since, for k ^ r, the 
function Fk{z), z £ G, is analytic in the interior of Ir U l[. and continuous on 
IrDl'j., we can replace in (I4.12p the path of integration Ir by with suitable 
orientation, i.e., for z £ G, 



-E 



(4.13) 



^ - z) 
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Observe that, by construction, /o is continuous on Zjv and Fk is continuous 
on In U for k ^ r. Thus, the function h in (14.13P is analytic in Gg and by 
Walsh's maximal convergence theorem there exist a sequence of polynomials 
{tn}^=i such that, 

\\h-tn\\L^(^G)^^^ '^eN, (4.14) 

where 1 < i? < ^. Since we can choose (i arbitrarily close to zi, ()4.14p is 
valid for any 1 < R < g. 

The function g in M.llh consists of sums of integrals of the type, 

9k{t) 
(t - z) 

where in view of (j4.4p and (j4.9p we have , for t £ Ij^, 

\9k{t)\ ^ |i-Tfcpffe+i. 

Hence, by using the result of Lemma 6 in [1], in conjunction with the re- 
mark following Theorem 1 of the same paper and the triangle inequality, we 
conclude that there exists a sequence of polynomials {qn}^=i satisfying 

y - qjLHG) ^ (4.15) 

where s := min{(2 — aii)^p''J_^_i : k = 1,2, . . . ,M}. This, combined with 
(j4.10p . (j4.14p and the triangle inequality, yields 

M ^ ^ 

\\fo-T.^k-(*- + 'lr.)\\mG)<ci-+C2 — , nGN. (4.16) 

k=l ^ 

Note that s = s*, if 7p^+i ^ N for the index k for which the minimum is 
attained in the definition of s. In the opposite case, where for the same 
index k, it holds 7p^+i G N, we get s = in (|4.16p by simply subtracting 
from g{z) and adding to h{z), in the right hand side of (j4.10p . the deriva- 

>fe + l^"^"'^ 

This observation and (jl.Sp imply that there exists a sequence of augmented 
polynomials {pn}, with pn G '^n^, such that, 

||i^(-,^o) -^«||l2{g) < C3;^ + C4-^, neN, (4.17) 

and the rest goes in similar lines as in the proof of Theorem 13. H except here 
we use the version of Andrievskii's lemma for functions with anti-derivatives 
in the space IP^^, given in [12^ Corollary 2.5]. These yield, 

1 r. 1 



five of the Cauchy integral on 1^, with density function a[,'^L i 

{z - rfc)^Pfc+i. 



Il/o - 7rn|lioo(G) ^ C5\/logn— + ceVlogn— , n > 2, (4.18) 
and ()4.7p follows with a different, but still less than g, R. □ 
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Remark 4.3. Note that < ii", as n — )• oo. Therefore from ()4.7p we 
recover the result of |12l Thm. 3.1]. However, Theorem 14.21 above gives, in 
addition, the precise dependence of the BKM/AB error on the pole singu- 
larities of /o for small values of n. We also note the lower estimate 

ll/o -7r„||ioo(G) > c;^, 71 e N, 
established in [12l Thm. 3.2]. 

4.3. BKM/AB with pole and corner singularities. We consider now 
the application of the BKM/AB with both pole and corner singular basis 
function of the form studied in Sections 13.21 and 14.21 Regarding poles we 
recall, in particular, our assumptions in Section 13.21 That is, the nearest 
singularities of /o in O are k poles, each one of order kj at zj, j = 1, 2, . . . k, 
where < |^>(z2)| < ••• < while the other singularities of /o 

occur at points z^+i, ^;k+2, • • •, where < |$(^ft+i)| < |$(2;k+2)| ^ • • • • 

Therefore, for the BKM/AB we consider the system defined by: 

(i) the K pole singular functions (j2.6p . with mj = 1, j = 1,2, . . . , k; 

(ii) the rM = Ylk=iPk corner singular functions of the form (j2.7p : 

(iii) and the n monomials (|2.8p . 
Accordingly, we let P^3 denote the space, 

K+rM+n 

{P ■ P{z) = ^ tjr]j{z), tj e C}, 
i=i 



1^3 

■ n 



and note that the associated augmented polynomial Kn{z,zo) is the best 
approximation to K{z, zq) in L^{G) out of the space P^^. 

The following result is a version of Andrievsii's lemma for functions with 
anti-derivatives in Pj^^. It will be used below, in the proof of the conclud- 
ing theorem of this section (in the transition from the L^(G)-norm to the 
L°°(G)-norm, where we establish the BKM/AB error in approximating /o 
by the augmented polynomials 7r„ derived from P^^. 

Lemma 4.1. Assume that T is piecewise analytic without cusps and let 
tk £ T , k = 1,2, . . . ,m. Also, let P G P„ and q be a fixed polynomial with no 
zeros on G. Assume further that for some constants an,k,j, k = 1,2, ... ,m, 
j = 1,2, . . . , Tk, the function 

P(z) ™ 

■' k=l j=l ^ 

where f ^(k) ^ {z) = {z—tk) j , with /3j > non-integer, satisfies: L{zq) = 
and \\L'\\l'^(^g) — ^ ■ Then, 

||i|lioo(G) <^^Myi^, (4.19) 
where C is a constant independent of n and o/ {{an,fc,j}yLi}fcLi- 
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Proof. The proof is based on Andrievskii's lemma for singular algebraic func- 
tions given in [121 Corollary 2.5] and relies on the results contained in [12^ 
§2]. The details of the derivation are as follows: 

First, we note that our assumption implies that F is a quasiconformal 
curve. Then, it is straightforward to verify that the results of Theorems 2.1 
and 2.2 (and hence the result of Corollary 2.2) in |12j hold true for functions 
of the form q'^{z)fj3^r{z), where fp^riz) := {z — r)^, with r € T and /? > 
non-integer. That is, 

inf -p||l2(G) - >2-l(0+i) ^ (4-20) 



n 

where avr (0 < a < 2) denotes the interior angle of F at r. 

With (|4.20p at hand it is, again, straightforward to verify consequentially 
that the results of Theorem 2.3, Corollaries 2.3 and 2.4, Lemma 2.3 and 
Corollary 2.5, of [12], hold true if we replace /^^ by f'^ ^. In particular. 
Corollary 2.5 of [12] applied to the function 

S{z) := / q^{z)L'{z)dz, 
where the path of integration [zq, z] is any rectifiable arc in G, gives that 



\\S\\loo(g) < ciVlogn ||5"||i2(G). 

(Note that S{zq) = and S'{z) = q^{z)L' {z).) Therefore, our hypothesis on 
||L'||^2(g) yields the inequality 



II^IIl-(g) ^ caVk^M. (4.21) 
On the other hand we have, 

L{z)= r q-\z)S'{z)dz = q-^(z)S{z) + 2 q-^{z)q'{z)S{z)dz, 

J Zq J Zq 

which implies 

and (j4T9]) follows from (fOT]) : cf. [H p. 122]. □ 

The concluding result of this section provides the theoretical justification 
for the use of the BKM/AB, with both corner and pole singularities. 

Let TTn denote the BKM /AB approximation to /o resulting from the space 
P^^. Then we have the following: 

Theorem 4.3. Assume that F is piecewise analytic without cusps and set 
Q := \<^{Zf,+i)\ and s* := min{(2 - 0^)7^,^^ : 1 < A; < M}. Then, 

_ I I 

ll/o - vrnllioo(G) ^ ciyOogn— + C2 — , n > 2, (4.22) 

for any R, 1 < R < g. 
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Proof. As in the proof of Theorem 13.21 we set Q{z) := Y[j=i{^ ~ Zj)^^~^^. 
The result (|4.22p will emerge by working as in the proof of Theorem 14.21 
The basic idea is to consider, in a bigger domain G, the anti-derivatives 
F and Gk of the functions QJq and Q-F^, respectively, in the place of the 
functions /q and Fj^.. The details of the derivation are as follows: 

We note that the function Q/q shares the same analytic properties with /g, 
apart from the fact that it has the singularities at the points Zj, j = !,...,«;, 
all removed. Therefore, the function 



can be extended analytically to a larger domain G than the one considered 
in Section 14.21 This larger domain G is obtained by choosing the point Ci 
close to the nearest pole z^+i of Q/q in but inside the level curve Lg, 
where now g := |^>(zk+i)|. The remaining part of the construction of G is 
exactly the same as in Section [4.21 

It follows therefore that (j4.23p is valid for z £ G, provided the arc of 
integration [zq,z] lies on G U dG \ {U^^-^r^} and is rectifiable. (This is 
always possible because dG is piecewise analytic.) Since the derivative of /o 
near can be obtained by termwise differentiation of the expansion (j4.4p . 
(cf. [3 p. 1448]) and since any power in the resulting expansion is bigger 
than — ^, we see that /q is integrable along any rectifiable arc in G with one 
endpoint at r^. Therefore, integration by parts gives, for z G G U dG, 



F{z) = Q{z)Uz)- Q'iOMOdC- Q'iOMOdC, (4.24) 



where we made use of the normalization of /o at zq. This shows that F is 
continuous on dG and analytic and on dG, except for the endpoints r^. By 
arguing as in (|4.8p we have, for z G G, 




(4.23) 






M 



(4.25) 



Similar properties to those of F apply to the anti-derivative 




(4.26) 



of QFl,k = l, 



...,M. That is, for z G G U 



Gk{z) = Q{z)Fk{z) - Q{zo)Fk{zo) 




(4.27) 
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and, for z £ G, 

M 

Next, by combining (|4.24p and (|4.27p we get. 



F{z)-dk-Gk{z) = Q{z)[Mz)-Fk{z)]- I Q'{C)[fo{z)-Fk{C)]dC, (4.29) 
where 



4 := Q{zo)Fk{zo) - r Q'(C)[/o(C) - FkiCM- 
This and ([33]) lead to 

\F{z) - 4 - Gk{z)\ ^\z- T^r^'+S z G 4. (4.30) 

By reasoning as in the proof of Theorem 14 . 2 1 we conclude, by using (|4.25p 
and (fOHjl that, for z E G, 

M 

Q{z)fo{z) - Q{z) Fl,{z) = g{z) + h{z), (4.31) 
k=l 

where the singular part 

of the splitting (j4.3ip can be approximated, eventually, by a sequence of 
polynomials {qn}'^=i at a polynomial rate, viz., 

\\S -qJiHG) ^ "-eN, (4.33) 



with s* := min{(2 — ak)^uj '■ 1 < A; < M} and the analytic part 
fc=l -^'fe 



" dk , If F(t) , 
dt + — - — ^-^dt 



(t - Zy 2TTi Ji^ {t - Z) 

}-Tri Z_' Z-^ /, (+ _ y\2 9„,- 



(4.34) 



can be approximated by a sequence of polynomials {tn}^^i at a geometric 
rate, viz., 

ll^-^nllLco(G) ^ nGN, (4.35) 
where 1 < < Hence using the triangle inequality we get 

M 

WQifo - E ^fc) - + 1n)\\LHG) <^l^+ ^2J^- (4-36) 
k=l 
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This implies 



M 

fe=l 



Q 



< 



(4.37) 



where d := 

niirij— ]^ 2,...,k{ I — ^j\ '• ^ ^ and <^ : — ^i* Thus, from 

p.5p we conclude there exists a sequence of augmented polynomials {i^n}) 
where p„ € P^^, such that, 

||K(-,Zo)-pn||L2(G) <C3;^ + C4-^, n G N, (4.38) 

Therefore, using the minimum property of the augmented kernel polynomi- 
als, we have 



||if(-,Z0) - i^n(-,2;o)||L2{G) < '^3;^ +C4-^, n G N, 

and this, in conjunction with the equivalence Lemma |2.H yields that 



ll/o - 7r^llL2(G) 



< 



1 



C5- 



71* 



+ C6 



1 



n G N. 



Since, 



vr„(2;) 



P(z) 



M Pk 



q(z) 



(4.39) 
(4.40) 
(4.41) 



k=l j=0 

where q{z) := 11^=1 (-^ ~ ^j)'^^ P{z) is a polynomial of degree n + ^, the 
rest goes as the concluding part of the proof of Theorem l4.2l except here we 
use the result of Lemma l4. II in the place of [121 Corollary 2.5]. □ 

5. NUMERICAL RESULTS 

Li this section we present numerical examples, that illustrate the con- 
vergence results predicted by the theory of Sections [3] and [H regarding the 
following four errors: 

E.,4K,G) := \\K{;Zo)-Kn{;Zo)\\L2^G), (5-1) 

En,oo{fo,G) := ll/o - ''rnllioo(g), (5.2) 

En,2{K,G) := \\K{;Zo)-Kn{;Zo)h^G), (5.3) 

En,oo{fo,G) := ll/o - 7fn|lx,oo(G)- (5-4) 

We do this by considering two different geometries: (a) lens-shaped do- 
mains; and (b) circular sectors. In both cases the normalized conformal map 
/o, and hence the kernel function K{-,zq), are known explicitly in terms of 
elementary functions. In addition, we present results illustrating the decay 
of the two sequences of points {Pnizo)}^=i and -Pn(-zo)^i of the Bergman 
polynomials. 
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5.1. Computational details. Let denote the set of linearly indepen- 
dent functions defined in (i2:6]) - (i2:8]) . For the apphcation of the BKM/AB 
(or BKM), we compute the associated orthonormal set {Pj} by using the 
Arnoldi variant of the Gram-Schmidt (GS) process studied in [20], rather 
than the conventional GS, which is based on the orthonormalization of the 
monomials {z^}, as it is suggested in and [16j . In the Arnoldi GS we 
construct first the polynomial part of the set {Pj} by orthonormalizing con- 
sequently the functions 1, zPq, zPi, . . . , zPn-i- Then, we orthonormalize the 
singular basis functions ()2.6p and ()2.7p . As it is shown in |20j , in this way we 
avoid the instability difficulties associated with the application of the con- 
ventional GS method. For a comprehensive report of experiments testifying 
the instability of the conventional GS in BKM and BKM/AB we refer to 
[Ml §5]. 

The GS process, requires the computation of inner products of the form 



{Vk,m) = / Vkiz)m{z)dA{z). (5.5) 
Jg 

For our purposes here, we compute these inner products by using Green's 
formula in order to transform the area integral into a line integral. For 
instance, when r]k = z^ , r]i = z\ we have 

In all cases considered below this leads to explicit formulas for the inner 
products (j5.5p . 

Regarding the computation of the errors (|5.ip - (j5.4p we note the following: 

(i) The two errors ||K(-, zo)-i^'„(-, zo)||l2(g) and \\K{-, zo)-Kni-, zo)\\l2(^g) 
are computed by using Parseval's identity, i.e., 

\\K{-,zo) - K„(-,zo)||^2(G) = K{zo,zo) - Kn{zo,zo), (5.7) 

and 

\\K{-,zo) - K„(-,zo)||^2(G) = K{zo,zo) - Kn{zo,zo). (5.8) 

(ii) Estimates for the two errors ||/o — T^nWi-x^Q^ and ||/o ~ ^Tn 1 1 ^oo^fj) are 
obtained by using the exact formula for /o and then sampling the 
differences /o — vr„ and /o — TTn on 100 uniformly distributed points 
on each analytic arc forming the boundary T. 

All results were obtained with Maple 11, using the systems facility for 
64-digit floating point arithmetic, on a pentium PC. 



5.2. BKM and BKM/AB approximation. 
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5.2.1. Lens-shaped domains. Let Ga,b denote the lens-shaped domain, whose 
boundary F consists of two circular arcs Ta and Tb that join together the 
points i and —i {Ta being to the left of Tb) and form angles a and b, respec- 
tively, with the hnear segment [— i, i]. (Thus, with the notation of Section [2?T] 
we have ai = 02 = ot, where a := {a + b)/7r.) Let /o denote the normahzed 
conformal map from Ga,b onto D{0,rQ), with /o(0) = and /o(0) = 1. By 
working as in [131 §4], it is easy to check that, ii a + b = /cvr/m, where 
/c, m G N, then /q is given by 



fo{z) = ro 



(-1)- 



Z+'l 



Z G Gah, 



-D-e-^'"- 



(5.9) 



where ro = {k/m)sm{ma/k). Also, 



K{z,0) 



and thus 



4m 



[iz-i)iz + i)]^ 



vr 



k { — i) k [z — i) fc 



' k [ij k [z + i) k 



2 ' 



(5.10) 



K{0,0) :-- 



1 



vr/c^ sin^(ma/A;) 

It is also easy to verify that the formulas (71)-(73) of [TJ] work as well 
for the exterior conformal map <^ : C\Ga,b — ^ ^ consider here. That is, 
w = $(z) is given by the composition of the following three transformations: 



az 

i(^) := ^m/{2m-k)^ 

wit) : 



^i{{m~k)-K / m+a) ' 



1 



arg^ G 



z ^ i 

-k-K jra^ {2m — k)TT/m], 

i((m—k)iT+ma)/{2m—k) 



(5.11) 

(5.12) 
(5.13) 



t-\a 

We consider separately the following three cases: 

(i) a = 1/2, with o = 7r/6 and b = 7r/3; 

(ii) a = 1/2, with a = 7r/4 and b = '7r/4; 

(iii) a = 2/13, with a = 7r/13 and b = 7r/13. 

Cases (i) and (ii): In the first two cases the conformal map /o is a rational 
function, and hence it has an analytic continuation across T into When 
a = 7r/6, then the two nearest singularities of /o in Vl are the two simple poles 
at zi = -\/3/3 and Z2 = \/3, where |<I>(2;i)| ^ 1.347 and \^{z2)\ ~ 2.532. 
Accordingly, in our experiments, we use the singular function \\-/{z — zi)]' . 
This cancels out the nearest singularity at zi. In the symmetric case, where 
a = 6 = 7r/4, we have 

and the only singularities of /o are the two simple poles at zi = — 1 and 
^2 = 1, where |$(zi)| = |$(z2)| = \/3. In this case, we use the singular 
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function \zj{z^ — -Zi)]', which takes care of both poles at z\ and 22- It 
follows from Remark 13. II that this cancels out all the singularities of /q. 
We recall from Theorems 13. II and 13.31 (and their proof) the four estimates, 



< En oo(/o, G) ^ "^y^ , (5.15) 

and 

— Tl 

^-^(^'^)^^(i^' ('-^^^ 

' ^^n.oo(/o,G)^!#f^. (5.17) 



Below, we present numerical results that illustrate the laws of the above 
errors and rates. In presenting the numerical results we use the following 
notation: 

• Q : This denotes the order of approximation (the base of n) in the 
errors (l53il) - (|5T711 . 

• Qn '■ This denotes the estimate of g, corresponding to n, and is deter- 
mined as follows: With denoting any of the two errors En^2{K, G) 
or En^2{K, G), we assume that 

Tl 

En ~ c— (5.18) 
and seek to estimate g by means of the formula, 

_L_ 

Sn=( " (5.19) 
\n-m En J 

(Here we take m = 4, or m = 5.) If En denotes either of the two 
errors £'n,oo(/oi G) or En^Qo{fo,G), then we assume that 

En - c!^, (5.20) 
and seek to estimate g by means of the formula, 

(5.21) 



1 



n -y/logn En 



yn-m y/log{n - m) En J 
with m = 4, or m = 5. 

gn '■ With En denoting either of the errors -E„,oo(/oi G) or -E'n,oo(/oi G), 
we also test the law 

En-c^, (5.22) 
thereby estimating g by means of 
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n 


BKM: 1.347 


BKM/AB: t 


2 w 2.532 




Qn 


EnAK, G) 


Qn 


5 


4.4e-01 




2.7e-02 




10 


1.3e-01 


1.47 


3.6e-04 


2.72 


15 


3.5e-02 


1.41 


4.1e-06 


2.65 


20 


8.9e-03 


1.39 


4.6e-08 


2.60 


25 


2.2e-03 


1.38 


4.9e-10 


2.59 


30 


5.4e-04 


1.37 


5.2e-12 


2.57 


35 


1.3e-04 


1.36 


5.4e-14 


2.57 



Table 5.1. BKM approximations to K: Lens-shaped, Case (i). 



n 


BKM: g 


« 1.347 


BKM/AB: 


2.532 


En,oo{fo, G) 


Qn 


Qn 


En,oo{fo, G) 


Qn 


Qn 


5 


2.5e-01 






1.4e-02 






10 


6.8e-02 


1.299 


1.54 


1.3e-04 


2.541 


3.03 


15 


1.6e-02 


1.331 


1.47 


1.3e-06 


2.528 


2.79 


20 


3.8e-03 


1.342 


1.43 


1.2e-08 


2.537 


2.71 


25 


8.5e-04 


1.346 


1.42 


l.le-10 


2.532 


2.67 


30 


1.9e-04 


1.347 


1.41 


l.le-12 


2.532 


2.64 


35 


4.3e-05 


1.347 


1.40 


l.le-14 


2.532 


2.62 



Table 5.2. BKM approximations to /q: Lens-shaped, Case (i). 



The presented results show clearly the advantage of the BKM/AB over 
the BKM. In addition, they indicate a close agreement between the theoret- 
ical and the computed order of approximation. In Tables 15.11 and 15.31 the 
results associated with the errors En^K^G) and En,2{K,G) indicate the 
convergence of gn to g. Regarding the errors En,ooifo,G) and En^ooifo^G), 
the results of the Tables 15.21 and 15.41 show that converges faster to g 
than gn- This suggest, at least for the geometry under consideration, a be- 
havior of the type ()5.22p for the errors En^oo{fo,G) and En,oo{fo,G). As 
it is predicted by Remark 13. 1^ in Case (ii) the two errors En^K^G) and 
En,oo{fo, G) vanish. This was testified in our experiments, in the sense that 
the computed errors En,2{K^G) and -E'„,oo(/oi G) were zero within machine 
precision, thus they are not quoted in Tables 15.31 and 15. 4[ 

Case (in): In this case the conformal map /o has a branch point singularity 
at each of the two corners ti = i and T2 = —i, and therefore Lehman's 
expansions (|2.4p are valid with 7^"'^^ = 'y^^ = 13/2 and 72^^ = ^!^^ = 1 -|- 



28 



M. LYTRIDES AND N. STYLIANOPOULOS 





BKM: Q ^ 1.732 


n 




Qn 


4 


2.7e-01 




8 


4.0e-02 


1.92 


12 


5.3e-03 


1.83 


16 


6.7e-04 


1.80 


20 


8.3e-05 


1.78 


24 


l.Oe-05 


1.78 


28 


1.2e-06 


1.77 


32 


1.4e-07 


1.76 


36 


1.7e-08 


1.75 



Table 5.3. BKM approximations to K: Lens-shaped, Case (ii). 



n 


BKM: g ^ 


1.732 




En,oo{fo, G) 




Qn 


4 


1.3e-01 






8 


1.6e-02 


1.685 


2.11 


12 


1.8e-03 


1.718 


1.95 


16 


2.0e-04 


1.729 


1.89 


20 


2.3e-05 


1.732 


1.83 


24 


2.5e-06 


1.732 


1.83 


28 


2.8e-07 


1.732 


1.81 


32 


3.1e-08 


1.732 


1.80 


36 


3.4e-09 


1.732 


1.79 



Table 5.4. BKM approximations to /q: Lens-shaped, Case (ii). 



l/a = 15/2. This gives (2 -a)/a = 12 and (2 - a)(l + I /a) = 180/13 = 
13.84 • • • . Furthermore, it follows from ()5.9p that the nearest singularities of 
/o in ri, are the two simple poles at zi = tan(7r/13) and Z2 = — tan(7r/13), 
where |$(2;i)| = |<I>(2;2)| ~ 1.119, and the next singularity occurs at a point 
Z3, where \^{z'i)\ w 2.055. 

Therefore, from Theorem 14.11 we have that, 

^n,2(K,G)^Cl^+C2-|^, (5.24) 

and 

I \ 

En,oo{fo, G) < csv^logn— ^ + ^47^7, (5.25) 

where 1 < R < \^{zi)\. In order to decide which singular functions to in- 
clude in the BKM/AB the following estimates, valid for n = 32, are relevant; 
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see also Theorem 14.31 



^ 8.7 X 10"^^ ^ 2.7 X 10-^ 



^ ~ 1.4 X 10"2\ w 1.0 X 10-^°. 



^(2-Q)(l+l/a) • ' |$(22) 

The estimates in the first hne indicate that for n = 32 (even for bigger 
values of n), the dominant term in the errors ()5.24p and ()5.25p is C2-^- As 
it is suggested by the estimate in the second line, we use in our BKM/AB 
approximations only the singular function [z/{z'^ — zf)]', which takes care 
of the two symmetric poles at zi and Z2, and we include no basis functions 
reflecting the corner singularities of /o on T. Then, from Theorem 14.31 we 
have for the resulting approximations that 

En,2{K,G) <ci^ + C2^, (5.26) 

En,oo{fo, G) ^ csTk^^I + ^4^' (5-27) 

where 1 < i? < |<1*(22)|- 

Below, we present numerical results that illustrate the rates in (15.240 - 
()5.27p . In presenting the numerical results we use the following notation: 

• Q : This denotes the order of approximation (the base of n) in the 
errors (iOill - ffOTl) . 

• Qn : This denotes the estimate of g, corresponding to n, and is 
determined as follows: With En denoting any of the four errors 
En,2{K,G), En,2{K,G), En,oo{fo,G) or En,oo{fo,G) we assume that 

En^c^, (5.28) 
and seek to estimate g by means of the formula, 

(5.29) 



1 

En 



The results quoted in Tables [531 and [5^ show the remarkable approxima- 
tion achieved by the BKM/AB by using as little as 32 monomials. Moreover, 
they highlight the significance of Theorem 14.31 &s it is compared to the es- 
timate (I1.13p . in the sense that they confirm fully the theoretical prediction 
that the two poles at zi and Z2 are the most serious singularities of /o for 
small values of n; see also Remark I 



5.2.2. Circular sector. Let Ga denote the symmetric circular sector of radius 
2 and opening angle an, < a < 2, at the origin, i.e., 

Ga ■= {z : l^l < 2, — Q7r/2 < argz < a7r/2}. 



30 



M. LYTRIDES AND N. STYLIANOPOULOS 



n 


BKM: Q « 


1.119 


BKM/AB: q ^ 2.055 


En,2{K, G) 


Qn 


EnAK, G) 


Qn 


4 


2.8819 




7.3e-02 




8 


2.3812 


1.049 


5.6e-03 


1.898 


12 


1.3864 


1.145 


3.9e-04 


1.934 


16 


0.9188 


1.108 


2.6e-05 


1.965 


20 


0.5961 


1.114 


1.7e-06 


1.974 


24 


0.3812 


1.118 


l.le-07 


1.982 


28 


0.2413 


1.121 


7.2e-09 


1.981 


32 


0.1538 


1.119 


4.6e-10 


1.992 



Table 5.5. BKM approximations to K: Lens-shaped, Case (iii). 



n 


BKM: Q ^ 


1.119 


BKM/AB: g 


PS 2.055 


-E'n,oo(/0) G) 


Bn 


En,oo{fo, G) 


Bn 


4 


0.8820 




1.2e-02 




8 


0.3817 


1.232 


6.1e-04 


2.101 


12 


0.2044 


1.170 


3.4e-05 


2.053 


16 


0.1180 


1.147 


2.0e-06 


2.029 


20 


0.0702 


1.139 


1.2e-07 


2.028 


24 


0.0424 


1.135 


7.0e-09 


2.028 


28 


0.0259 


1.131 


4.1e-10 


2.028 


32 


0.0160 


1.128 


2.5e-ll 


2.028 



Table 5.6. BKM approximations to /q: Lens-shaped, Case (iii). 



Let /o denote the normahzed conformal map from Ga onto D(0, tq), with 
/o(l) = and /o(l) = 1- For each value of the parameter a the confor- 
mal map fo{z) can be computed by means of the transformations (see |12^ 
p. 532]): 



where 



This gives 



''o 



iz 



2a(4i/° - 1) 



iz^/a _ 2l/a 



41/a + 1 



and d 



t - d 
td-l 



2I/C 



2l/a 



2a(4V" - 1) 
41/" + 1 



and K{1,1) 



1 / 41/a + 1 



TT V2a(4V« - 1) 



(5.30) 



(5.31) 



(5.32) 
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The normalized exterior map <l> : C\Gq, — A is given, as can be easily 
verified, by the composition of the following three transformations: 

m ■■= e'\ arg^ G (-^/2,3^/2], (5.34) 
1 - e^'^/H 

w{t) := (5.35) 

We consider separately the following two cases: 

(i) a = 1 (half-disk); 

(ii) a = 3/2 (three-quarter disk). 

Case (i): When a = 1, then the domain Ga is the half-disk 

Gi = {z: \z\ <2,d\z> 0}. 

In this case the conformal map /o has an analytic continuation across T into 
Q. The nearest singularities of /o in 0, are the two simple poles at zi = —1 
and Z2 = 4, where |$(2;i)| ~ 1.452 and |$(z2)| ~ 2.212. Accordingly, in our 
experiments we use the singular function [l/{z — zi)]', which cancels out the 
nearest pole at zi. This case is similar to the lens-shaped domain with a = 
1/2. Hence, the errors£^,2(^^), ^n,oo( /o, G) , En^ {K, G) and ^„,oo(/o, G) 
satisfy respectively (j5.14p . (j5.15p . (j5.16p and (j5.17p . Our purpose here, is 
to illustrate that the error bounds in (j5.14p - (|5.17p reflect the actual errors. 
We do so by computing estimates to Qn and of g by using (j5.18p - (j5.23p . 

In Table [5?71 the results associated with the errors En,2{K, G) and En,2iK, G) 
indicate clearly the convergence of Qn to g. Regarding the errors En^ooifo, G) 
and En^oo{fo,G), the results of Table ESI show that £»* converges faster to 
g than This suggest a behavior of the type (I5.22p for En,ooifo,G) and 
En,oo{fo, G). In both tables the numbers confirm the remarkable advantage 
of the BKM/AB over the BKM. 

Case (ii): In this case /o has a branch point singularity at the point ti = 
with 

oo 

/o(^) = /o(0)+5^a,-^^/", ai^O, 
i=i 

valid for z close to 0. The nearest singularity of /o in is a simple pole 
at zi = 4, where |$(zi)| ~ 2.04. For the application of BKM, Theorem 14.11 
gives that 

^",2(i^,G)^ci-i^ + C2-^, (5.36) 

and 

\ \ 

En.ooifo, G) ^ C3 Vlogn— — - + C4 — , (5.37) 

where 1 < R < \^{zi)\. 

Since l/|<I>(2;i)|^'' ~ 3.3 x 10"^^, and in view of Theorem 14.31 we include 
in our basis only singular functions that reflect the branch point singularity 
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n 


BKM: 1.452 


BKM/AB: < 


?w 2.212 


E„ o{K, G) 


On 


En,2{K, G) 


Qn 


5 


7.1e-02 




2.8e-02 


- 


10 


1.6e-02 


1.55 


7.2e-04 


2.39 


15 


2.8e-03 


1.54 


1.7e-05 


2.29 


20 


5.1e-04 


1.49 


3.6e-07 


2.29 


25 


8.7e-05 


1.49 


7.6e-09 


2.26 


30 


1.5e-05 


1.48 


1.6e-10 


2.24 


35 


2.4e-06 


1.47 


3.2e-12 


2.24 


40 


4.0e-07 


1.47 


6.4e-14 


2.24 


45 


6.6e-08 


1.47 


1.3e-15 


2.23 


50 


l.le-08 


1.46 


2.6e-17 


2.23 



Table 5.7. BKM approximations to K: Half-disk. 



n 


BKM: Q 


^ 1.452 


BKM/AB: 


2.212 


En,oo{fo, G) 


Qn 


Qn 


En,oo{fo, G) 


Qn 


Qn 


5 


l.le-01 






4.0e-02 






10 


2.2e-02 


1.401 


1.64 


7.7e-04 


2.20 


2.62 


15 


3.4e-03 


1.446 


1.60 


1.5e-05 


2.20 


2.42 


20 


5.3e-04 


1.450 


1.55 


2.8e-07 


2.22 


2.37 


25 


8.3e-05 


1.452 


1.53 


5.2e-09 


2.22 


2.34 


30 


1.3e-05 


1.452 


1.51 


9.8e-ll 


2.21 


2.31 


35 


2.0e-06 


1.452 


1.51 


1.8e-12 


2.22 


2.30 


40 


3.1e-07 


1.452 


1.50 


3.5e-14 


2.20 


2.27 


45 


4.8e-08 


1.452 


1.49 


6.6e-16 


2.21 


2.27 


50 


7.4e-09 


1.452 


1.49 


1.2e-17 


2.22 


2.27 



Table 5.8. BKM approximations to /q: Half-disk. 



of /o at T\. More precisely, in order to keep the contribution of both somxes 
of error balanced, we choose to use the first 15 singular function of the form 
where j /a ^ N. This gives s* = 23/3 in Theorem 14.31 ^-^id hence the 
fohowing estimates for the errors in the resulting BKM/AB approximations, 

^«.2(i^,G)<ci-^ + C2-^, (5.38) 

and 

En,oo{fo,G) ^ Csv^log 71-23^3 + — , (5.39) 

where 1 < R < \^{zi)\. 
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Below, we present numerical results that illustrate the rates in (I5.38p - 
(j5.39p . where we use the following notation: 

• a : This denotes the exponent of 1/n in the errors (I5.38p - (j5.39p . 

• an : This denotes the estimate of cr corresponding to n, and is deter- 
mined as follows: With En denoting any of the two errors En^2{K, G), 
En,2 {K, G) , we assume that 

En ~ (5.40) 
and seek to estimate a by means of the formula 

a„ = logf^Vlogf^y (5-41) 



^ En J \n 

If En denotes either of the two errors -E'n,oo(/o) G) or -En,oo(/o) G), 
then we assume that 

En ~ c^bi^^, (5.42) 
and seek to estimate a by means of the formula 

an = ^" ; " • (5.43) 

In addition, we check a behavior of the form (j5.28p for the errors En,2{K, G) 
and En,ooifoi G), by computing Qn as in (j5.29p . with 5 in the place of 4. 

Our purposes here is to show that the change of the dominant term in 
both (|5.38p and (|5.39p can actually be detected in the computed errors. 
This is indeed the case in the results quoted in Table 15.91 More precisely, 
the results associated with the errors En,2{K,G) and En,oc{fo,G) indicate 
the convergence of Qn to g for values of n up to 50 and the convergence of 
(Tn to a for values larger than 50. Furthermore, the results show that the 
two constants ci and C2 in (j5.38p and C3 and C4 in (j5.39p are, respectively, 
of the same magnitude. 

5.3. Rates of decrease of the Bergman polynomials. First, we present 
results illustrating the rate of decrease of the sequence {P„(l)} for the cir- 
cular sector considered in Section 15.2.21 with a = 2/5. In this case, the 
nearest singularities of /o in Q are the two simple poles at the symmetric 
points zi = e^*'^/^, Z2 = e~^"^/^, where |$(2;i)| = |^(-Z2)| ~ 1.145. From the 
proof of Corollary 13.11 and (|5.36p we have that 

\Pn{zo)\ < \\K{;Zo)-Kn{;Zo)\\mG)<Cl^^^^^ + C2^, (5.44) 

where 1 < R < \^{zi)\, Accordingly, we check to detect the decay in the 
following two forms: 

\Pn{l)\^c^ and \Pn{l)\^c^, (5.45) 
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n 


BKM/AB 


: a r 


« 7.67 g « 2.04 








Qn 






Qn 


20 


7.2e-05 


- 


_ 


8.2e-05 


- 


_ 


25 


1.6e-05 


6.74 


1.35 


1.5e-05 


7.62 


1.40 


30 


z.ye-06 


y.oi 


1.41 


z.De-Uo 


y. /u 


1.42 


35 


2.2e-07 


16.84 


1.67 


1.8e-07 


17.37 


1.71 


40 


l.Oe-08 


22.81 


1.84 


7.7e-09 


23.44 


1.86 


45 


4.1e-10 


27.41 


1.90 


2.8e-10 


28.17 


1.94 


50 


1.3e-ll 


32.83 


1.99 


l.Oe-11 


31.25 


1.95 


55 


7.5e-12 


5.84 


1.12 


5.3e-12 


7.05 


1.14 


60 


2.6e-12 


11.84 


1.23 


2.0e-12 


11.20 


1.22 


65 


1.3e-12 


8.68 


1.15 


9.9e-13 


8.73 


1.15 


70 


7.4e-13 


7.50 


1.12 


5.9e-13 


7.04 


1.11 


75 


4.4e-13 


7.58 


1.11 


3.5e-13 


7.57 


1.11 


80 


2.7e-13 


7.66 


1.10 


2.1e-13 


7.62 


1.10 



Table 5.9. BKM approximations to /o and K: 3/4-disk. 



with Q = \^[zi)\ and, in view of the remark made in [121 pp. 530-531], 
a = {2 — a) /a + 1/2. We do so, by estimating g and cr, respectively, by 



means of the formulas 



|^n.^l0(l) 
\Pn{l)\ 



(5.46) 



-n = log(%i#Vlogf^). (5.47) 



and 

rr„ = lop- I J— , ^ , , , , 

|p„,(i)| J' 

The results listed in Table 15.101 show clearly the transition from one dom- 
inant term to the other in ()5.44p for values of n around 50. 

We end, by presenting results that illustrate the rate of decrease of the 
augmented sequence {P„(l)}, for the circular sector considered in Section 
15.2.21 where now we consider the two cases a = 3/4 and a = 4/5. When 
a = 3/4, then the nearest singularities of /o in are the two simple poles 
at the symmetric points zi = ^"^'^1'^ and = e~^*'^/^, where |^>(zi)| = 
|^('2^2)| ^ 1.349. When a = 4/5, then the nearest singularities of /o in 
f2 are the two simple poles at the symmetric points zi = e,^™!'^ and Z2 = 
g-4i7r/5^ where |<I'(-2;i)| = |'&(-Z2)| ~ 1.372. In both cases, we construct the 
sequence {P„(z)} by augmenting the monomial basis functions with the 
singular function which reflects the branch point singularity of /o at 

Ti = 0, and we seek to detect the decay of the sequence {P„(l)} in the form 
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35 



n 


a = 


4.5 


1.145 


\Pn{l)\ 




Qn 


10 


2.6e-02 


- 


- 


20 


1.2e-03 


4.51 


1.37 


30 


7.6e-06 


12.38 


1.65 


40 


1.7e-06 


5.14 


1.16 


50 


4.0e-07 


6.57 


1.16 


60 


1.8e-07 


4.35 


1.08 


70 


9.1e-08 


4.49 


1.07 


80 


5.0e-08 


4.50 


1.06 


90 


2.9e-08 


4.50 


1.05 


100 


1.8e-08 


4.50 


1.05 



Table 5.10. Rate of decrease of |Pn(l)|: Circular sector, a = 2/5. 





a « 3.833 


n 


\Pn{l)\ 




10 


2.8e-03 




20 


7.2e-05 


5.30 


30 


l.le-05 


4.60 


40 


3.2e-06 


4.36 


50 


1.3e-06 


3.86 


60 


6.6e-07 


3.89 


70 


3.6e-07 


3.89 


80 


2.2e-07 


3.88 


90 


1.4e-07 


3.88 


100 


9.1e-08 


3.87 



Table 5.11. Rate of decrease of |P„(1)|: Circular sector, a = 3/4. 



where, in view of Theorem 14.31 and [T2l pp. 530-531], a = 2(2 — a)/a + 1/2. 
As above, we estimate a by means of the formula 

-n = logf%#Vlogf^y (5.48) 



,P„(1)| J \n-lQ 

The results listed in Tables 15.111 and 15.121 indicate clearly the conver- 
gence of (Tn to the predicted value of u, indicating that the argument in [12 
pp. 530-531] applies also to the case of the augmented Bergman polynomials 
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a = 3.5 


n 




^ Tl 


10 


2.3e-03 




20 


2.7e-04 


3.10 


30 


2.2e-05 


6.17 


40 


8.7e-06 


3.20 


50 


3.8e-06 


3.72 


60 


2.0e-06 


3.66 


70 


l.le-06 


3.63 


80 


6.9e-07 


3.61 


90 


4.5e-07 


3.59 


100 


3.1e-07 


3.58 



Table 5.12. Rate of decrease of |P„(1)|: Circular sector, a = 4/5. 
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